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Abstract 

Is it possible to trap a quantum particle in an open geometry? In 
this work we deal with the boundary value problem of the station- 
ary Schroedinger (or Helmholtz) equation within a waveguide with 
straight segments and a rectangular bending. The problem can be 
reduced to a one dimensional matrix Schroedinger equation using 
two descriptions: oblique modes and conformal coordinates. We use 
a corner-corrected WKB formalism to find the energies of the one- 
dimensional problem. It is shown that the presence of bound states is 
an effect due to the boundary alone, with no classical counterpart for 
this geometry. The conformal description proves to be simpler, as the 
coupling of transversal modes is not essential in this case. 
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1 Introduction 

The use of symmetry methods to understand bound states has an ancient 
origin and has given many of its fruits in atomic and subatomic physics - 
The work of Marcos Moshinsky is an excellent example of this pp. Over 
the years, however, other approaches to understand more complicated sys- 
tems have given up symmetry in both classical and quantum mechanics [2J, 



[3J. Here we present one of the simplest ways of destroying continuous sym- 
metries with purely quantum-mechanical consequences: A corner. We will 
show that a sharply bent waveguide supports a finite number of bound states 
whose nature has no classical counterpart. Due to the geometry, the tradi- 
tional methods of integrability cannot be used to obtain the corresponding 
solutions. 

This type of systems has been studied in the past, either with mathemat- 
ical tools [H El El E] proving the existence of bound states or by employing 
numerical methods for computing the corresponding energies and eigenfunc- 
tions [51 El HD] • The approach in our considerations will be completely ana- 
lytical, albeit we use approximations to perform computations. We want to 
show two different approaches to the problem in order to exhibit an appar- 
ent complexity of the description given in other studies, mostly numerical [8]. 
First we shall deal with a convenient basis of states which shows the binding 
as a consequence of the presence of corners: A sharp obstacle diffracts waves 
of an arbitrary wavelength. Some limitations of this method, such as the 
coupling of an infinite number of modes, will be pointed out. Finally, the use 
of conformal coordinates will be introduced as a natural description to find 
all the bound states for arbitrary bending angles. 

We would like to emphasize that interesting applications may come in 
the form of bent wave guides carved in crystals, as suggested in [TTj, us- 
ing methods connected to the evolution of Airy profiles - a recent example 
can be found in [12] . The landscape looks promising, considering the many 
applications of controlling, bending and focusing waves. 

2 The existence of bound states 

2.1 The oblique modes and the effective potential 

We show how to obtain a one dimensional Schroedinger equation with an 
effective potential along the longitudinal direction of the guide. One of the 
challenges for describing our system is to find a suitable set of states. A 
natural way to introduce box eigenmodes in our system is by using oblique 
modes, as shown in the figure. 

Let us take units such that h 2 /2m = 1 with m the mass of the particle 
and write the stationary Schroedinger equation in the form 
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Figure 1: Description of our system. The distance between corners is denoted 
by I. The slope a is related to the bending angle a. We show an oblique 
mode with n = 2 



VjLV + £ty = (!) 



For the configuration in figure 2.1 with Dirichlet conditions, we express the 
solution in the oblique basis as 

^( x , y) = Y, sin ( ~r(y ~ a \ x \ + ) ( 2 ) 

n=l ^ 1 ' 

and ij}(x,y) ^ for a|x| — I < y < a\x\. The matrix elements of the 
Schroedinger operator in this basis can be obtained by using the expan- 
sion Q and integrating over the y variable. Let us denote the vector wave 
function by ^f(x). Using the following definition for matrix elements 

TflTi 

N =5 n M = (1 - (-) m+n ) — (3) 

and employing a gauge transformation given by 

= exp ^- 2C y^M^J (4) 
N(x) = expl p M ) Nex P I — i M j ( 5 ) 
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we reduce the Schroedinger equation to 



d 2 



71 



- (-) (a 2 + l)N(a;) 2 - ~ M 2 + E 



2a \ 2 
/ ) 



dx 2 V I 

where we see an effective potential of the form 



&(x) = 



(6) 



V 



eff 



[a 2 + l)N(x) 2 + 



2a\ 2 2 



(7) 



This matrix couples, in principle, all components of The origin of M 
and N 2 is related to the momentum and the energy of a particle in a box. 
The x dependence of the potential V e ff comes completely from the non- 
commutability of N and M and is non-differentiable at the corner due to \x\. 
Another important feature is that a change of units x' = x/l,E' = El 2 in 
equation (|6f leads to a scale-free equation, thus proving that the energy E 
must scale as l~ 2 . Therefore, the presence of bound states in this problem is 
independent of the scale. This is in contrast with other bent systems with 
smooth boundaries (5j EJ U\ ■ 

2.2 Solutions from corner-corrected WKB 

In previous works by Wheeler, Bestle and Schleich [H], the WKB approx- 
imation was corrected in the presence of potentials with corners. As a re- 
sult, the Schroedinger equation with a potential could be approximated by 
a WKB equation with a source term. By neglecting completely the effect 
of the potential except at the non-differentiable point, the WKB equation 
with a source could be interpreted as a free Schroedinger equation with a 5 
potential. In our case, this would cause a binding effect depending on the 
energy-dependent amplitude of the 5. We shall apply this idea in order to 
cope with V e ff obtained in the last section, with the additional feature that 
V e ff is a matrix potential. 

We start with the Schroedinger equation 



d 2 u(x) 

dx 2 



(x)u(x 



(8) 
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where A = (E — V e ff) is now a matrix. The corner-corrected WKB 
formalism given in [H] can be easily generalized to matrix potentials. A 
WKB scattering equation (such as eq (13) in [13]) can be obtained for (Tsj) . 
The resulting equation, in turn, can be treated by keeping first order terms 
in a and considering a negligible potential except at x — 0. This gives the 
isolated effect of the corner in the form 



an 
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2 ,M] E-^N 2 ) 6(x). (9) 



The energies for bound states in this problem are obtained by i) imposing 
the appropriate boundary conditions u(x) = at x = ±oo, ii) substituting 
the corresponding solutions in the condition for the jump in the derivative 
(due to 5), namely 



du, du, an 



2 



dx 0+ dx I 3 
The solutions are found to be 




\N\M] [E--N 2 u . (10) 



7T 



-1/4 



u(x) = Af I — 



E 



cxp 



-N 2 



E u 



(11) 



with a normalization factor. The vector u is determined by the condition 



coming from the jump in the derivative (10), namely 



an 



[N 2 , M] + 2 



u = 



(12) 



This is a linear equation for u and the wave function can be obtained by 



choosing uq as a member of the null space of the operator in (12). The scaled 



energies E = (l/n) 2 E are obtained from (12) by using the determinant 
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det 



7T 



[N 2 , M] 



N — E 



3/2 







(13) 



The 2x2 case can be solved easily by inserting in (13) the diagonal matrix 
N 2 



T2 _ ft'/* 



diag 



ground state is 



E^ 



^3/2 



. The resulting energy for the 




This energy is below the threshold (rr/l) 2 . The level E\ can also be obtained; 
it is above 4(7r/Z) 2 . Therefore, in this 2x2 approximation, the only state 
which decays exponentially with the distance is the ground state. The wave 
function can be found straightforwardly once Uq is obtained. We show two 
cases in figures |2]|3j 

The solutions thus obtained show that a bound state develops due to 
the corner alone, producing a 5 potential and coupling the oblique modes - 
such coupling is essential, since a 1 x 1 approximation makes the interaction 
vanish. Yet, we may ask: where are the excited states below threshold due 
to an increasing value of a? In the following, we analyze this problem using 
conformal coordinates as a way to obtain a suitable basis. Such a basis 
cannot be obtained by means of perturbation theory in a using the solutions 
of this section. 



3 The conformal map 

Our aim now is to show that a non-trivial boundary value problem such as 
the bent waveguide can be mapped smoothly to a straight guide with an 
effective potential (equivalently, a problem with position-dependent mass). 
The simplest way to describe this process is by using complex variables. 
Denote z = x + iy, £ = u + iv and consider a transformation such that 
£ = F(z). We impose the following conditions on F: 

1) The image of the bent waveguide is an infinite straight strip described 
by < u < F(c). 
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Figure 3: Ground state wavefunction for a = 0.2. The peak becomes more 
pronounced as a increases. For large values of a, the description is no longer 
valid. 

2) The Laplacian is transformed into a quadratic form without cross- 
terms (The matrix M will not appear). 

3) The transformation is smooth except at the points O (upper corner) 
and c (lower corner) in figure |4j 

The condition 1) allows the use of a complete set of functions {sin(nu/F(c))} 
for the Dirichlet problem and {cos((n + l/2)u/F(c))} for the Neumann prob- 
lem, given that the contour lines meet the boundaries orthogonally. The 
condition 2) precludes the appearance of terms proportional to d u : The in- 
tegral Jq F ^ sm(nu/F(c))d u sm(mu/ F(c)) would produce unbounded matrix 
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Figure 4: Contour lines of the coordinates u, v. We denote the lower corner 
by c and place the upper corner at the origin O, see the text. 

elements due to well-known properties of the momentum operator for a par- 
ticle in a one-dimensional box. Requirement 3) ensures that any suspicious 
behavior of the transformation is due to the corners. 

3.1 Explicit transformation 

The requirements described above are fulfilled by choosing F analytic for all 
the points in the interior of the waveguide. Orthogonality is reached by the 
Cauchy-Riemann conditions. We introduce the operators 

d z = d x -id y , d z *=d x + id y (15) 

= d u - id v , <%* = d u + id v (16) 

obtaining thus = d z d z *, V\ v = d^*. Using the analytic function 
£ = F(z), the Laplacian operator is transformed as <9g<9£» = \F'(z)\~ 2 d z d z *. 
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Using the Cauchy-Riemann conditions we have |F'(z)| 2 = d(u,v)/d(x,y), 
i.e. the Jacobian of the transformation as a function of z. Our boundary 
value problem with the Schroedinger equation becomes 



<%9 e + | — \ 2 Ej^(u,v) = 0, ue[0,F(c)}, veR (17) 

We give an explicit example of our conformal map in the following. Consider 
the Schwarz-Christoffel (SC) transformation [T5] of the upper half-plane onto 
the semi-infinite strip. By using a composition of a SC map for a trigon with 
internal angles 7r/2 together with a SC map for a trigon with internal angles 
air and (1 — a)ir, we arrive at the bijective transformation between one half of 
the bent waveguide and a semi-infinite straight strip. The expression results 
in 



rs'm 2 (z/2) rz 

£ = F(z) = / dx(l - x)- a x~ a = / d0(tan(0/2)) 9 

Jo Jo 

= /3(sin 2 (z/2),^,i± g ) (18) 

where we have used the definition of the incomplete beta function [16] and 
we have defined q = 1 — 2a as the complementary angle of 2a in 7r radians. 
Therefore < q < 1. In passing we note that F(c) = sec(7rg/2). The 
Jacobian is given by ||| 2 = |cot^/2| 29 and one can approximate it in terms 
of £ (or u,v) for some cases of interest: 



dz f/(?)l£l 2q ^ q+1 ^ near the origin 

l^l 2 = 9(q)\F(c) - ei 2 ^ 1 "^ near F(c) (19) 
[ 1 + 4g cos^e - !"! away from the corners 

with f(q) = ((1 + q)/2)- 2q K q+1 \ g(q) = (2(1 - q)) 2 ^^. It is important 
to note that as one moves away from the corner in the direction of the 
arms, the Jacobian (19) becomes a hyperbolic cotangent of a real variable 
and approches to unity exponentially fast. Therefore, interactions vanish for 
points far from the corners. In connection with the existence of other maps, 
we stress that any other transformation constructed in this way will have 
the same behavior near the corners: The fractional power law of the map 
is completely determined by the angle of the guides and therefore by our 
physical system. 
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3.2 Approximate solutions in conformal coordinates 



One can try many methods of solution for the transformed boundary value 
problem given above. Here we give an example based on the fact that our 



careful construction allows certain approximations. Let us reduce (17) to a 
one-dimensional problem by defining e n = —n 2 /F(c) 2 and 



T _ . ,,dz. 2 . 2F(c) r*F(c) ( nu\ , ( mu\ dz 2 



The set J2F(c)/n sm(nu/F(c)) is obviously orthonormal. After multiply- 



ing (17) by our complete set of functions and performing the corresponding 



integrals over u, we find that (17) is equivalent to 



(Snmdl - EV nm (v)) 1p m (v) = 6 n 1p n (v) 



(21) 



Near the origin, our potential (20) becomes 



It is justified to approximate the potential by its diagonal entries V n (v), since 
V nm (v) decreases as we move away from the diagonal for fixed v. This can 
be seen by noting that for n ^ m the integrand can take negative values and 
oscillates rapidly asn-m increases. We have 



(-d 2 v - EV n (v)) $ n (v) = e n ^ n {v). (23) 
We include some plots of V n (v) and its derivative for different angles, figures 

EE 

The potentials shown in the figure have zero derivative at v — (except for 
q — 1) but, apart from that, their behavior is that of a shallow potential well 
whose derivative suffers an abrupt change in a small region near the corner. 
In the plot showing the slopes, we see that there is almost a jump in V^(v) 
around v — 0. The formalism developed by Wheeler, Schleich and Bestle [H] 
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Potential V(v) for different values of q 

V in units of E/e 



-0.6 

^^^rS^^S^ -0.8 




VV^^^r--^ 1 ' 




\ \^1> 




q=0 


-1.6 


q=l/2 





Figure 5: A set of potentials for n — 1, q increases from to 1/2 in 10 steps. 
The transversal mode energy is shown at —1. 



Slope V'(v) for different values of q 

V in natural units 



q=3/4 
q=l/2 

q=l/4 





Figure 6: The slope of the potentials as a function of v. Near the origin we 
see an abrupt change (a discontinuity appears in the limit q — 1) 
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in order to deal with 'kinks' should be helpful here again. Although the only 
case in which the derivative is strictly discontinuous corresponds to q — 1, 
the strong variation of V^(v) occurs in a region which looks effectively like 
a corner for long wavelengths (we expect this approximation to work at low 
energies). As the simplest approximation, let us consider the corner-corrected 
WKB formula for quantization of bound states, namely 



S = V2 £ dv\Jt n + EV n (v), Classical orbit action (24) 



AS = 2(k + |)tt - 4arctan(^), Phase shift (25) 



A = E (K(0+) - K(0-)) (2(e n + EV n {0))y 3/2 , Jump in the slope (26) 
Remarkably, one can obtain good results by neglecting the classical orbit 



S = and keeping the corner effect alone. Using (24), the approximation 
implies 



A 
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1, k even (27) 



finding thereby that all even states obey this relation, but the excited states 
are not supported by such a weak potential. The excitations come from the 



solutions of (27) by varying n and not k. This results in a transversal node 



structure, instead of longitudinal. With these considerations, the solution to 



our problem (27) requires V n (0) and V^(0) as the only input. Such quantities 



can be estimated by means of (22) for m = n. As mentioned before, we 
replace the jump in the derivative V^(0+) — V^(0— ) by twice the maximum 
of Vn(v) around the origin (see the figure). We denote such quantity as V^ ax . 



Squaring (27) we obtain a cubic relation for E : 



(n 2 cos 2 (7rg/2) - EV n (0)f = ^E\V^ ax f (28) 
One can estimate V n (0) and V^ ax as 



12 



Energy levels as a function of the angle 

n=l n=2 




Figure 7: Energ 1GS clS db function of the angle. There is a bound state for all 
q > as a consequence of the corner. An unlimited number of bound states 
appear as q increases. 



V n (0) ~ (cos(7rg/2)) 29/(y+1) , 



V^ ax ^\q{ C o^q/2)f^ +l \ (29) 



With this, the energies in (28) can be expressed in terms of n and q using 



Cardano's formula. We present the results graphically in figure [7j 

The agreement of the estimated solutions with previous numerical results 
[TOl [T7] is qualitatively good, considering that the action has been completely 
neglected in favor of the phase shift. Also, the off-diagonal elements of the 



matrix potential (22) have been ignored. The trend of the curves as functions 



of the bending angle is comparable with the numerical results obtained in 
[TO] for several modes: the energies decrease and more states appear below 
threshold as q h-> 1. The node structure of our wavefunctions is also compat- 
ible with the solutions obtained in [TU] both numerically and experimentally. 
Remarkably, the fractional power behavior of the Jacobian alone is respon- 
sible for the structure of the spectrum of this system. We include a plot of 
the resulting wave density for the ground state of an L-shaped guide. 

In what touches the corrections, we expect the following: The action 



appears to be important as q increases (deeper wells). Therefore (28) receives 
corrections which increase the energy. As a consequence, we expect a shift 
of the critical angles towards the value q = 1. Precisely at this value of q 
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Figure 8: Contourplot of the probability density for q = 0.5 and n = 1 
(ground state energy E ~ 0.56). The peak of the distribution is located near 
the upper corner, as expected. This is achieved with only one transversal 
mode in conformal coordinates 

we also expect a non-zero lower bound for the energy E m i n ~ 1/4, as can be 
seen by identifying such an extreme case with a guide of double thickness. 

4 Discussion 

In the first part we have shown that the existence of bound states could be 
understood as a corner effect and the coupling of transversal modes. We could 
learn that diffraction was essential to the final result, but the coupling of the 
modes gave an extra complication. The second part, however, shows that 
an appropriate set of coordinates may keep the effect of the corner without 
the coupling of such modes. This makes contact with recent dicussions [18] 
in which perturbation theory seems doomed in the presence of unbounded 
curvature - our example fits perfectly in such class of problems. One can 
follow either approach to solve the problem, but its apparent complexity 
- e.g. the many modes used in [8] for the L-shape - dissapears using the 
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second method. At the end we can conclude that a simple way of breaking 
a continuous symmetry such as a corner still allows for a simple description 
as long as the important information is used: diffraction and conformality. 

The next task is to find the correction to the energy levels due to the 
classical orbits entering through the action. But we can say that for now 
our goal has been reached. As a final comment, the problem of bound states 
in multiple junctions -e.g. X-shaped and *-shaped geometries - lends itself 
for this kind of treatment by using Neumann boundary conditions. Due to 
the properties of our conformal coordinate system, this can be achieved by 
replacing sm(nu/F(c)) by cos((n + \)u/ F(c)) in our discussion. 
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